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RESUMEN 

Las estrellas de neutrones poseen los campos magneticos mas intensos que conocemos. Aquf se pretende 
dar una discusion pedagogica de parte de la ffsica relacionada. Las estrellas de neutrones existen gracias al 
principio de exclusion de Pauli, en dos sentidos: 1) Este hace diffcil comprimir muchas partfculas en un espacio 
pequeno, permitiendo un estado de equilibrio mecanico en presencia de una fuerza de gravedad muy intensa. 
2) La ocupacion de estados cuanticos de baja energfa por parte de protones y electrones impide el decaimiento 
beta de neutrones poco energeticos. Un corolario de esto es que necesariamente hay partfculas cargadas en 
la estrella, permitiendo el flujo de corrientes electricas. Como las partfculas son degeneradas, colisionan muy 
poco, haciendo posible la cxistcncia prolongada de campos magneticos intensos. Estos se revelan en pulsarcs 
y son la mas probable fuente de energfa para la alta luminosidad en rayos X y gamma en "magnctarcs" . 
Discuto brevemente el posible origen de estos campos magneticos, asf como algunas consideraciones ffsicas que 
restringen a sus configuraciones de equilibrio. 

ABSTRACT 

Neutron stars have the strongest magnetic fields known anywhere in the Universe. In this review, I intend to 
give a pedagogical discussion of some of the related physics. Neutron stars exist because of Pauli's exclusion 
principle, in two senses: 1) It makes it difficult to squeeze particles too close together, in this way allowing 
a mechanical equilibrium state in the presence of extremely strong gravity. 2) The occupation of low-energy 
proton and electron states makes it impossible for low-energy neutrons to beta decay. A corollary of the second 
statement is that charged particles are necessarily present inside a neutron star, allowing currents to flow. Since 
these particles are degenerate, they collide very little, and therefore make it possible for the star to support 
strong, organized magnetic fields over long times. These show themselves in pulsars and are the most likely 
energy source for the high X-ray and gamma-ray luminosity "magnetars" . I briefly discuss the possible origin 
of this field and some physical constraints on its equilibrium configurations. 

Key Words: dense matter — (magnetohydrodynamics:) MHD — magnetic fields — stars: magnetic fields — stars: 
neutron — stars: white dwarfs 



1. INTRODUCTION 

This text aims at giving a pedagogical intro- 
duction to the physics of neutron stars and their 
magnetic fields, starting from undergraduate physics 
(Quantum Mechanics and Electromagnetism) and 
leading up to some current research questions such as 
the origin and equilibrium configuration of the mag- 
netic field and clarifying some misconceptions ap- 
pearing in the recent literature. All equations used 
here can be derived "on the back of an envelope" , 
based only on undergraduate physics, and students 
are strongly encouraged to do so. 

In the philosophy of keeping the treatment sim- 
ple and aiming at a good understanding of basic 
questions rather than introducing all the complica- 
tions that might conceivably arise, I consider an ex- 



tremely conservative model of neutron stars that in- 
cludes neutrons, protons, and electrons, all treated 
as degenerate and mostly non- interacting fermions. 
This neglects many extremely interesting (but un- 
certain) issues such as neutron superfluidity, pro- 
ton superconductivity, or quark deconfincmcnt and 
condensation. For more comprehensive treatments 
of neutron star physics, also addressing many of 
t he "exotic" issues, I sugg est the classic book of 



Shapiro & Tcukolsky ( 1983) and the recent volu mes 
of iGlendennind (|200Cf) and lrlaensel et all (|2007l) . as 
well as other works mentioned in specific places of 
the text. For an inspiring popular history of compact 
sta r research, writ ten by one of its main participants, 



sec 



Thornd (jl994) . 



2. DEGENERATE FERMIONS 
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The lowest energy state of a system of fermions is 
the "Fermi sea" or momentum-space "Fermi sphere" , 
in which N fermions occupy the least energetic or- 
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bitals available to them. If the fermions are con- 
fined to a real-space volume V and otherwise non- 
interacting, they will agglomerate in momentum 
space in a sphere around p = with radius pf 
(the "Fermi momentum"), so the volume occupied 
in 6-dimcnsional phase space (x,y, z,p Xl p y ,p z ) is 
V 6 = V x (4np F /3). It can be shown (e. g., by 
solving the Schrodinger or Dirac equation in a box) 
that there are h~ 3 single-particle orbitals per unit 
phase-space volume, where h = 2itTi is Planck's con- 
stant. For spin-1/2 fermions such as electrons, pro- 
tons, or neutrons, Pauli's exclusion principle allows 
two particles to be put in each orbital (with spin 
projections s z = ±1/2 on an arbitrary axis), thus 
the total occupied phase-space volume must also be 
equal to Ve = Nh 3 /2. Equating both expressions for 
1/6, one finds the Fermi momentum as a function of 
the number density of fermions, n = N/V, 



PF 



^TT 2 ™) 1 / 3 , 



(1) 



which holds regardless of how relativistic the parti- 
cles are. 

Global properties of the fermion system, such as 
its total energy E, can be calculated as sums over all 
particles, i. e., integrals over momentum space, 



d 3 p e(p), 



(2) 



p\<PF 



where one should generally use the relativis- 
tic energy-momentum relation e(p) = [(mc 2 ) 2 + 
(cp) 2 ] 1 / 2 , where m is the mass of the fermions and 
c is the speed of light. These integrals can b e done 



analytically (see Shapiro fc Teukolskv 1983); how- 



ever, the calculation greatly simplifies in both the 
non-relativistic limit (p <C mc), in which e(p) « 
mc + p 2 /(2m), and in the extreme relativistic limit 
(p ^> mc), in which e(p) ~ cp. In these limits, the 
pressure P = —{dE/dV)M takes a polytropic form 
Pan 7 , respectively 
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and 
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(3) 



(4) 



(For rough estimates, it is useful to remember that, 
aside from a factor 0.4 — 0.25, the pressure is just n 
times the kinetic part of the Fermi energy.) 

3. WHITE DWARFS AND NEUTRON STARS 

Since the temperature is essentially zero (com- 
pared to the Fermi energies) , the structure of degen- 
erate stars is determined by only two of the standard 



four equations of stellar structure, whose Newtonian 
(weak-gravity) form is 



dP GM{r)p dM(r) 



dr 



dr 



Anr 2 p, (5) 



which can be combined to yield an order-of- 
magnitude expression for the pressure required to 
sustain a star of given total mass M and character- 
istic mass density p against its own gravity, 



P ~ GAf 2 / 3 // 3 . 



((») 



The ratio p/n is generally constant, taking simi- 
lar values: 

• p/n w Am u /Z w 2m u in white dwarfs, in which 
most of the mass density is provided by ions of 
mass Am u (where m u is the atomic mass unit) 
and charge Ze, while the degeneracy pressure is 
due to the electrons, and 

• p/n ~ m n w m u in neutron stars, where neu- 
trons (of mass m n ) dominate both the pressure 
and the mass density. 

In the low-density, non-relativistic limit of 
eq. ([3]), the latter can be combined with eq. (J6j) to 
estimate an equilibrium radius for the star, 



R 



Gm 



-AT 



(7) 



Due to the widely different fermion masses m (m„ « 
1839 m e ), the radii for white dwarfs and neutron 
stars are very different (~ 10 4 km for the former 
and ~ 10 km for the latter, at a similar, solar mass 
M ~ A/q). Eq. ((7|) shows that, as the mass increases, 
the stars become smaller and denser and thus their 
matter more relativistic. In the limit of high den- 
sity, the hydrostatic pressure of eq. ([6]) needs to be 
provided by extremely relativistic particles (eq. [4]), 
yielding a maximum mass 



( mm 

V m u 



M, 



0! 



(8) 

where mpi = (hc/G)? is the "Planck mass", 
a natural mass scale for quantum gravity. For 
white dwarfs, this maximum mass is the well-known 
"Chandrasckhar limit" Aichandra ~ 1-4 Af , beyond 
which a degenerate stellar core collapses to form a 
neutron star or black hole. Note that the formation 
of a neutron star involves the release of the bind- 
ing energy ~ GM 2 /R ~ 0.2 Mc 2 as neutrinos (and 
perhaps gravitational waves), so a neutron star mass 
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can be lower than the Chandrasekhar limit by a cor- 
responding amount. In fact, a neutron star mass as 
low as 1.25 Mq was measur ed in the recently discov- 
ered double pulsar system (jLvne et al 1 l2004h . 

White dwarfs are well described by Newtonian 
gravity, and the main inter-particle forces are elec- 
trostatic, for which the interaction energies are much 
smaller than the kinetic energies, so the physics de- 
scribed above suffices to obtain an accurate descrip- 
tion of the structure of a white dwarf. In neu- 
tron stars, such estimates are generally not accurate, 
since at super-nuclear densities the baryons (neu- 
trons and protons) interact strongly with each other 
(interaction energies comparable to the kinetic ener- 
gies), and the effects of General Relativity become 
important (see below). Instructions for students to 
construct more realistic neutron star models num eri- 
cally have been giv en bv lSilbar fc Reddvl (120041 ) and 



Sagert et al 



( 2006f ). Eq. ||SJ) suggests a maximum 
mass for neutron stars at most a few times larger 
than Mchandra- This is likely to be correct, but its 
precise value is not known, mainly due to the uncer- 
tain strong interactions among the neutrons. 

Given the basic structural parameters of the stars 
(here taken as M = Mq and R = 10 km), one 
can obtain several other interesting quantities. The 
escape speed v csc = (2GM/R)? = 0.5 c confirms 
that neutron star gravity is "strong" and general- 
relativistic effects can be important . Their minimum 
allowed rotation period (at which the centrifugal 
force on the equator equals the gravitational force) 
P min = 2n[R 3 /(GM)]i = 0.5 ms, is much shorter 
than for any other kind of stars (including white 
dwarfs, for which P m i n ~ 20 s), a crucial argument 
in identifying pulsars as neutron stars. 

Finally, one may estimate a safe upper bound 
on a typical magnetic field B in a neutron star by 
requiring that the magnetic energy be lower than the 
gravitational binding energy (so the Lorentz force 
does not exceed the gravitational force jl, 



B 2 AnR 3 GM 2 

< 



yielding 
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B„ 



R 



10 18 G. 



(9) 



(10) 



4. WHY NEUTRON STARS? 



The name "neutron stars" incorrectly suggests 
stars composed exclusively of neutrons. However, 
additional particles inside these stars play a crucial 



2 As customary in theoretical astrophysics, I use Gaussian 
cgs units, in which the unit of magnetic field is the gauss, 

1 _ 3 

G = erg 2 cm 2 . 



role. A neutron (n) in vacuum decays by the weak 
interaction process n — > p+e + i? e (beta decay) into a 
proton (p), an electron (e), and an electron antincu- 
trino (v e ), with a half-life close to 15 minutes. This is 
impeded in very dense matter by the Pauli exclusion 
principle: If all the low-energy proton and electron 
states are already occupied, only suffiently energetic 
neutrons can decay. On the other hand, if many pro- 
tons and electrons are present, some of these will be 
energetic enough to combine into neutrons by inverse 
beta decay, p + e — > n + v e , where v e stands for an 
electron neutrino. In a neutron star, the neutrons 
and protons will be confined by gravity, the elec- 
trons by the electrostatic pot ential of the protons 
(e. g., Reisenegger et al. 20061 ). while neutrinos and 
antineutrinos are unbound and esc ape, contribut- 
ing t o the cooling of the star (e. g., lYakovlev et al. 
20011 ). Direct and inverse beta_ decays are in bal- 



ance if the chemical potential^] of neutrons (/J, n ), 
protons (n p ), and electrons (/x e ) satisfy the relation 
Mn = L l p + Me 7 which forces the coexistence of a 
small fraction (few percent, but density-dependent) 
of charged part icles with a much larger nu mber of 
neutrons (e. g., IShapiro fe Teukolsky Il983h . Addi- 
tional particles (both charged and uncharged) can 
appear by other weak decay processes at densities 
higher than typical nuclear densities. 

In addition to stabilizing the neutrons, the 
charged particles play two important roles regard- 
ing the magnetic fields and their evolution. Being 
charged, these particles can generate electrical cur- 
rents, which support potentially very strong mag- 
netic fields. In addition, since the proton fraction 
Y = rip/n depends on density (rii stands for the 
number density of particle species i = n,p,e, and 
n = n n + n p is the total baryon density), neu- 
tron star matter is inhomogeneous, stabilizing it 
with respect to convective overturn (Pethicklll992l 
Reisenegger fc Goldreich 1992 ; Reisenegger 2001al ). 



As discussed below, this is likely to have an impor- 
tant stabilizing effect on magnetic field configura- 
tions. 

5. FARADAY'S LAW AND ASTROPHYSICAL 
MAGNETIC FIELDS 

Long-lived magnetic fields are ubiquitous in the 
Universe, and neutron stars are no exception. These 
nearly static magnetic fields must have currents act- 
ing as sources, according to Ampere's law, 



VxB = — j. 

c 



(11) 



3 At zero temperature, these chemical potentials reduce to 
the respective Fermi energies. 
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Currents imply charges of one sign (i. e., electrons) 
moving with respect to those of the opposite sign 
(protons or other ions) . On time scales much shorter 
than the observed lifetimes of the fields, these par- 
ticles suffer Coulomb collisions, which would damp 
their relative motion, if it were purely due to inertia. 

Thus, in order for the magnetic fields to survive, 
a much larger "inertia" is required, which is provided 
by Faraday's induction law: Any change of a mag- 
netic field induces an electric field, 



dB_ 
~dt 



-cV x E. 



(12) 



The electric field keeps the current going in spite of 
the frequent collisions, as described by Ohm's law, 



l = E, 
a 



(13) 



where the conductivity a is inversely proportional 
to the collision rate. Combining eqs. (fTTj) . (fT2|) . 
and (|13p . one obtains a diffusion equation for the 
magnetic field, 



v = -Vx| —V 
dt \4ir<T 



(14) 



suggesting an Ohmic (or resistive) decay time ts ~ 
4iraR 2 /c 2 , where R is a characteristic length scale, 
such as a stellar radius. In astrophysical plasmas, the 
length scales are enormously larger than in labora- 
tory conditions, which allows astrophysical magnetic 
fields to be so long-lived. In the particular case of 
neutron stars, the Pauli principle also makes it hard 
for particles to be scattered into different quantum 
states, therefore enhancin g the conductivity. Based 
on this, Bavm et al. ( 19691 ) showed that neutron star 
magnetic fields might live longer than a Hubble time. 

The only "loophole" that allows for a significant 
evolution of the magnetic field is to allow for a ve- 
locity field v e of the electrons in the reference frame 
of interest (i. e., the center-of-mass frame of a star). 
This velocity adds a magnetic force term (l/c)v e x B 
on the right-hand side of eq. (|13[) . and correspond- 
ingly a term V x (v e x B) on the right-hand side of 
cq. (|14p . which can be interpreted as an advection 
of the magnetic field lines by the motion of the elec- 
trons. In most astrophysical contexts, this motion is 
shared by all the other particles, correspond ing to an 
ideal magneto-hydrodynamic motion (e. g., Kulsrudi 
2005h . H owever, in neutron stars, other variants are 
possible ( Goldreich fc Reisenegger Il992l ). such as a 
motion of all the charged particles with r espect to the 
neutrons ("ambipolar diffusion", e. g.. iHovos et al 



20081 ). or by only the most mobile charg e carrier 
("Hall drift"; e. g. . iReisenegger et"aL 2007 ). each of 
which has quite distinctive properties. 

6. MAGNETIC FIELDS IN NEUTRON STARS 

Many neutron stars are detected as pulsars, 
whose regular pulsations in the radio, X-ray, and/or 
optical bands are produced by a strong magnetic 
field turning around at the stellar rotation period 
P. These periods slowly increase in time, i. e., 
the neutron stars lose rotational energy, probably 
through magnetic coupling with their surroundings. 
Modelling this coupling as electromagnetic radiation 
from a dipole rotating in vacuum, oriented orthog- 
onally to the rotation axis, one can infer the sur- 
face magnetic field strength B oc \J PP, where P 
is the time-derivative of th e rotation period (e. g., 
Shapiro fc Teukolskvl Il983l ). Inferred fields range 
from 10 8 G in millisecond pulsars up to 10 15 G in soft 
gamma-ray repeaters (SGRs), the latter being the 
strongest magnetic fields known in the Universe, but 
still 3 orders of magnitude weaker than the strongest 
that might conceivably be present in neutron stars 
according to eq. (fTU|) . In this (dynamical) sense, neu- 
tron star magnetic fields are quite weak, as they are 
in all other stars known so far. In spite of this, the 
magnetic field may be the main agent breaking the 
axial symmetry of the mass distribution in a rotat- 
ing a neutron star, in this way producing preces- 
sion (Wasserman 2003), as appears to be observed in 
some pulsars, and gravitational waves, which might 
quickly reduce the rotation rate of newborn neutron 
stars (Cutler 2002). 

In radio pulsars, the rotational energy loss can 
account for the whole observed energy output (rela- 
tivistic particles and electromagnetic radiation). For 
the strongly magnetized, but slowly rotating SGRs 
and anomalous X-ray pulsars (AXPs), however, the 
observed X-ray luminosity is much larger than the 
rotational energy loss rate, so an additional source 
of energy is required, t he most likely being the deca y 
of their magnetic field ([Thompson fc Duncanl[l996h . 
This would make these objects be the only known 
magnetically powered stars, or "magnctars" . An in- 
teresting way of probing the strong magnetic fields 
inside these objects appears to be the quasi-periodic 
oscillations recently detected following two large 
flares of SGRs and interpreted as crustal shearing 
modes coupled to Alfyen w aves travelling through 
the stellar core ([Levinl 120071 ). 

In very old neutron stars, such as millisecond 
pulsars and low-mass X-ray binaries, the magnetic 
field is < 10 9 G, weaker than in young neutron stars, 
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such as radio pulsars and high-mass X-ray binaries 
(~ 10 n ~ 14 G), suggesting that the magnetic field 
strength decays with time, perhaps induced by ac- 
cretion of matter from the binary companion (e. g., 
Payne fc Melatosl2007 and references therein). Mag- 



netic field decay within the population of single ra- 



(Ostriker & Gunnlll969t Naravan & Ostriked Il990h 


but disputed by others ( 


Bhattacharva et al.l 


Il992t 

seem 


Faucher-Giguere & Kaspi 


20061), and does nol 



to be well established. 

7. ORIGIN OF THE MAGNETIC FIELD 

A natural hypothesis to explain the origin of 
the strong fields observed in neutron stars is the 
compression of the magnetic fl ux already pre sent 
in the progenitor stars. It led IWoltierl (Il964h to 
predict field strengths of 10 14 " 16 G, before any 
neutron stars had b een identified ob s ervationally. 
Many authors (e. g., Rudermanl 1972 :_ Reisenegger 
2001b : Ferrario fc Wickramasing hel l2005al lbl. l2006h 
have pointed out that the distribution of magnetic 
fluxes is very similar in magnetic A and B stars, 
white dwarfs, and neutron stars, in this way provid- 
ing support for the hypothesis of the magnetic fluxes 
being generated on or even before the main-sequence 
stage and then inherit ed by the compact remnants. 

On the other hand, Thompson fc Duncan (1993) 
suggested that newborn neutron stars are likely to 
combine vigorous convection and differential rota- 
tion, making a dynamo process likely to operate in 
them. They predicted fields up to 10 15 ~ 16 G in 
neutron stars with few-millisecond initial periods, 
and suggested that such fields could explain much of 
the phenomenology associate d with SGRs and AXPs 
( Duncan fc Thompson 19921 : Thompson fc Duncan 
19951 ll996h T some of which were later confirmed 
to spin down at a rate consistent with a stron 
dipole field Q0 14 - 15 G; 



Woods et al.lll999h . 



Kouveliotou et al. 199 



Of course, the two processes are not mutually 
exclusive. A strong field might be present in the 
collapsing star, but later be deformed and perhaps 
amplified by some combination of convecti on, dif- 



ferential rotation, and magnetic instabilities (jTavle 



1973t ISpruitll2002f ) . The relative importance of these 
ingredients depends on the initial field strength and 
rotation rate of the star. For both mechanisms, the 
field and its supporting currents are not likely to be 
confined to the solid crust of the star, but distributed 
in most of the stellar interior, which is mostly a fluid 
mixture of neutrons, protons, electrons, and other, 
more exotic particles. 



8. PERSISTENT, ORDERED FIELD 
STRUCTURES 

The magnetic fields of neutron stars, like those of 
upper main sequence stars and white dwarfs, appear 
to be ordered (with a roughly dipolar external con- 
figuration) and persistent (for much longer than a 
solar cycle, perhaps for the entire existence of these 
stars). As mentioned above, their magnetic flux dis- 
tributions arc similar, which also implies that their 
ratios of magnetic to gravitational energy (or mag- 
netic stress to fluid pressure) are similarly small in all 
of them. The Lorentz force is much smaller than the 
pressure gradient and the gravitational force that are 
dominant in establishing the hydrostatic equilibrium 
in the star. Therefore, the hydro-magnetic equilib- 
rium state can be considered as a small perturba- 
tion to an unmagnetized, "background" hydrostatic 
equilibrium (denoted by a subscript "0"), in which 
the (conceptual) introduction of the magnetic field 
forces the fluid to displace from its "initial" posi- 
tion, r —> r + £(r). This causes small perturbations 
of density, which are customarily described in two 
complementary ways: 

• Eulerian perturbations, which compare the den- 
sity at the same point in space, before and after 
the perturbation, 



8p(r) = p(r) -p (r), 



(15) 



• Lagrangian perturbations, the change in the 
same fluid element before and after being dis- 
placed, 

Ap(f) = p(r + ar\)-p (r). (16) 
These two descriptions are related by 

Ap-Sp = ^V Po = ^ r , (17) 

and exactly the analogous relations hold for the Eu- 
lerian and Lagrangian perturbations of the pressure, 
SP and AP. These perturbations must satisfy the 
force balance condition 



3 x B 



VSP - Sp Vipo = 0, 



(18) 



where ipo( r ) is the background gravitational poten- 
tial (assumed to be unperturbed, in the so-called 
"Cowling approximation"). 

Another important, shared property of these 
stars is that a large part (if not all) of their interior is 
stably stratified (i. e., it resists convective overturn). 
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In the case of upper main-sequence envelopes and 
white dwarfs, this is because of a radially increas- 
ing entropy; in the case of neutron stars, because 
of a radial dependence in the fraction of protons, 
electrons, and possibly other particles. This has the 
consequence that the adiabatic sound speed, relating 
the Lagrangian pressure and density perturbations 
of a given fluid element (that conserve entropy and 
composition), c 2 s = AP/Ap, is larger than the back- 
ground derivative dPo/dpo = (dPo/dr)/(dpo/dr), in 
which entropy or composition are changing (e. g., 
Reisenegger fe Goldreichl Il992h . The Lagran gian 
perturbations can be directly related to the diver- 
gence of the fluid displacement field £(r) causing 
them, 

AP 

— = Ap = -p V • £■ (19) 

From all the above, one obtains that the Eulerian 
perturbations Sp and SP are linearly independent 
combinations of V • £ and £ r , 



5p = -p V • £ - ^£r, 
dr 

6P=-p c 2 s V-Z-^ r , 



(20) 
(21) 



and can therefore be regarded as independent vari- 
ables. 

Thus, a given, sufficiently weak magnetic field 
B(r) corresponds to an equilibrium configuration if 
two independent scalar functions SP and Sp can be 
found that satisfy eq. (fl"8|) . Since the latter is a 3- 
component vector equation, this will not generally 
be possible. Therefore, it imposes a condition on B 
that can be written as 



V x [(V x B) x B] = 0, 



(22) 



a single, scalar condition on the magnetic field. 
This condition is much less restrictive than the of- 
ten imposed force-free cond ition, (V x B) x B = 
( Broderick fc Naravan|[200 8), relevant for the oppo- 
site limit of dynamically dominant fiel ds, or even the 
cond ition V x [(V x B) x B/p] = (jHaskell et al ' 



20071 ). required only for barotropic fluids with a 
unique pressure-density relation p(P), which are not 
stably stratified. 

The assumption of axial symmetry, with 



B = Vx [a(r, 0)V<j>] + /3(r, 6)\7<p, 



(23) 



where V0 = </>/(rsin#), considerably simplifies the 
problem of constructing equilibria, since in this case 
cq. (|18p implies that the Lorentz force can have no 



azimuthal component, (j x B)^ = 0, therefore sur- 
faces of constant (3 coincide with those of constant 
a (i. e., (3 = (3[a(r, 9)]). These "allowed" magnetic 
field configurations produce two independent force 
components in the meridional plane, which can gen- 
erally be cancelled by an appropriate choice of the 
functions SP and Sp, so no additional conditions are 
imposed on B(r) to correspond to an equilibrium. 

Of course, in order to be viable, a certain mag- 
netic field configuration must correspond to a stable 
equilibrium, which is much more difficult to charac- 
terize and still largely an open problem. Analytic at- 
tempts to search for stable magnetic field configura- 
tions have failed, only yielding the general result that 
both purely toroidal fields B = (3(r, 6*)V0 and purely 
poloidal fields B = V x \a(r, fl)V < / >1 are unstable 



(|Tavlerl 119731 ; iFlowers fc Rudermanl 119771) . and the 
speculation that linked toroidal and poloidal fields 
might st abilize each other, yielding a st able equi- 
librium ( PrendergastJ 1956 ; Wright 19731). Recent 
MHD simulations dBraithwaite fc Spruitll2004 [2005 



Braithwaite fc Nordlundl I2006T ) have shown initially 
complex, "random" magnetic fields to evolve on an 
Alfven-likc timescale into a roughly axisymmetric, 
linked poloidal-toroidal configuration that persisted 
for a resistive timescale and thus might be a good 
approximation to the field structures in upper-main 
sequence, white dwarf, and neutron stars. 

Once a stable, ideal-MHD equilibrium mag- 
netic field has been established, it will sur- 
vive for many Alfven times, but not for- 
ever, since there are several dissipative processes 
by which it could evolve on long time scales 
dGoldreich fc Reiseneggerl 1992 ; Reisenegger 2007 



uvjQiareicn <y; n eisoncggor ryyz; jneisenegger i/uu/[ 
Reisenegger et all 2007t Hovos et al. 2008 ). possibly 
matching the times on which magneta r fields appear 
to decay ( Thompson fc Duncan 19961 ). 



9. CONCLUSIONS 

Neutron stars are fascinating objects with ex- 
treme properties, which include the strongest mag- 
netic fields in the Universe. Nevertheless, these share 
properties with those of other stars, among these, 
that they are weak in the sense of producing only 
small disturbances to the structure of the respective 
stars. Some progress has been made in understand- 
ing possible magnetic field configurations and their 
evolution, but there is still much left to do. 
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